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In this study, using g-generalized bit cumulants (g is the nonextensivity parameter 
^ , of the recently introduced Tsallis statistics), we investigate the asymmetric unimodal 

maps xt+i = 1 ~ a|a::t|^' (i = 1,2 correspond to > and Xt < Q respectively, Zi > 1, 
< a < 2, t — 0,1,2,...). The study of the g-generalized second cumulant of these maps 
allows us to determine, for the first time in the literature, the dependence of the infiexion 
' parameter pairs (01,02) to the nonextensivity parameter q. This behaviour is found to 

be very similar to that of the logistic-like maps {zi = Z2 = z) reported recently by Costa 
et al. [Phys. Rev. E 56 (1997) 245]. 
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c/j ' In recent years, sensitivity to initial conditions of nonlinear dynamical systems has started to be 
studied with an increasing interest. Among them, dissipative systems of low-dimensional maps 

' self-organized criticality gj, symbolic sequences |^ and conservative systems of long-ranged many-body 

^ Hamiltonians low-dimensional maps |^ can be enumerated. We focus here on the one-dimensional 

' , dissipative maps. As it has already been well studied in the literature, to investigate the sensitivity to 

^ initial conditions of any kind of one-dimensional maps on the onset to chaos, it is possible to introduce 

R ■ X , ^x(t) , , 



(where Aa;(0) and Ax(t) are discrepancies of the initial conditions at times and t) which, if the Lyapunov 
. exponent Ai 7^ 0, satisfies the differential equation d^^/dt = Ai^, thus 

5^: m^e^^' , (2) 



whereas, if Ai = 0, it satisfies d^/dt — A^^^, thus 

at) = [i + {i-q)Xgt]^ (qen) (3) 



which recovers the standard case (namely Eq.(2)) for q = 1. Here, q is the nonextensivity parameter 
of recently introduced Tsallis statistics and for g ^ 1 case, it is evident that Eq.(3) yields a power- 
I law sensitivity to initial conditions. At the onset of chaos, ^(t) presents strong fluctuations reflecting 
(-H ' the fractal- like structure of the critical attractor and Eq. (3) delimits the power-law growth of the upper 

O ' bounds of £,{t). These upper bounds (^ (x i~) allow us to calculate the values of q for the map under 
consideration. This method of findingnumerical values of q has been successfully used for the logistic 
map [0 , a family of logistic-like maps , the circle map |^ and a family of circular- like maps [Q . Beside 
this method, another one has been developed by Lyra and Tsallis |^ by looking at the geometrical aspects 
of dynamical attractors at the treshold to chaos. Using the multifractal singularity spectrum /(a) |pT[|, 
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they proposed the scaling relation 



(4, 



*e-mail: tirnakli@sci.ege.edu.tr 



1 



where otmax (ftmin) IS the most rarefied (concentrated) region of the multifractal singularity spectrum 
of the attractor. This relation presents the second method of calculating the q values once the scaling 
properties of the dynamical attractor are known. It has already been shown that for all above-mentioned 
one-dimensional dissipative maps, the values of q index calculated within these two different methods are 
the same (within a good precision) for any given value of the inflexion parameter z of the used map. In 
addition to this, another important information which comes from all these works is the determination 
of the behaviour of the q index as a function of the inflexion parameter z. As it can be seen from Fig. 2 
of 1^ and from the Table of |Q, it seems that when z ^ oo , q approaches 1. 

The purpose of this paper is to show numerically whether such kind of behaviour is verified for asym- 
metric unimodal maps (AUMs) 

_ J 1 - a\xtY^ if xt > , . 

^ \ l~a\xtY-^ if < 

where Zi^2 > 1, < a < 2 and i = 0, 1, 2, .... Unfortunately, up to now, for the AUMs, neither the first nor 
the second method can yield satisfactory results for the prediction of q values of any inflexion parameter 
pair (zi, Z2) [0. The problem of the first method might be originated from the prediction of the critical 
Oc values at the chaos treshold with enough precision, whereas the problem of the second method might 
be related to the numerical procedure used to estimate the /(a) curve, whose most rarefied region [amax] 
is usually poorly sampled. Another fact which yields this failure might be the nonuniversal behaviour 
of AUMs reported in [|l3|,^ (for example, it is shown that AUMs fail to exhibit the metric universality 
of Feigenbaum). Since the values of q are not available for AUMs, clearly it is not possible to see the 
behaviour of g as a function of (zi, Z2) pairs by using the above-mentioned two methods. In this paper, 
in order to see this behaviour, without finding the precise values of q for (21,^2) pairs, we use another 
technique based on the very recent generalization of bit cumulants for chaoti c syst ems within Tsallis 



statistics |l5|,|l^. To explain this technique, let us recall the main results of [y_5 16|. The generalized 
second cumulant (or heat capacity) is given by 
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{p"'-') - {P"?] (6) 



where p is the natural invariant density ||l7| ]. As 5 —> 1, we have 



Cl'> ^l^{\npY)-{\npr (7) 

which is equivalent to the standard definition of the second cumulant |T7|JT^ ]. It is shown that 

there is a kind of scaling between c'^'^ and Cj^"* which is evident from a €2'^ vs c'^'^ plot, where it is 
easily seen that most of the points can be fitted to a straight line. The slope of this line gives the scaling 
factor between c'^'^ and C^'^ and also provides us a tool to fulfill our aim in this paper. In Fig.l, we 

illustrate this slope plotting Cj^'' vs c'^^ curve for a representative (zi, Z2) pair. It is evident that most 
of the data points fall onto a straight line which yields a well-defined slope (in fact, we encountered that 
a few points deviate from this straight line especially when the values of the inflexion parameters start 
to be very different from the standard case (namely, zi = Z2), but since such data points are very few, 
we prefer to calculate the slope using a linear regression -with an error of ±0.01 for each estimation of 
the slope- to data points which fall onto this straight line). As it is evident from Fig. 2, when q ^ 1, 
naturally Cj"^^ C2^^ and thus the slope also tends to unity. Therefore, now we have another parameter 
(which behaves exactly the same as q) from where we can check the behaviour of the q index as a 
function of (^1,^2) pairs without knowing the exact values of q\ Calculating the above-mentioned slope 
for various values of (zi,Z2) pairs, we are able to estimate the behaviour of q index as a function of 
inflexion parameters without finding any q value. Fig. 3 represents this behaviour for (2, 22)7(2^1, 2), (3, Z2) 
and {zi, 3) cases. It seems from the figure that as Z2 — zi ±00, the above-mentioned slope (and thus, q 
index) approaches to unity] This tendency is exactly the same as that observed for a family of logistic-like 
maps and a family of circular- like maps Q, namely, as z — > cxd, q approaches 1. 

Summing up, for the first time (to the best of our knowledge), we manage to study the asymmetric 
unimodal maps in such a way that it became possible to relate it with the nonextensivity parameter q 
of Tsallis statistics. We were able to show that the dependence of the inflexion parameter pairs (zi, Z2) 
of these maps to the q index is very similar to that observed for other onc-dimcnsional maps reported so 
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far. Although the determination of the exact values of q index for (zi, Z2) pairs is still lacking, we hope 
that present work would be considered as a first attempt on this line and will accelarate other studies 
which are, no doubt, highly welcome. 
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Figure Captions 



Figure 1 : The scaling between the standard second cumulant C2 and the generalized second cumulant 
c'^^ for a representative value of {zi, Z2) pairs. 



Figure 2 : The behaviour of the slope as a function of q index for a representative value of (21,^2) 
pairs. 



Figure 3 : The behaviour of the slope as a function of Z2 — zi values for a family of four different (zi, Z2) 
pairs. 
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